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ABSTRACT 

Context. The Spitzer GLIMPSE and MIPSGAL surveys have revealed a wealth of details of the Galactic plane in the infrared (IR). 
We use these surveys to study the energetics and dust properties of the Eagle Nebula (Ml 6), one of the best known SFR. 
Aims. We present MIPSGAL observations of M16 at 24 and 70 yum and combine them with previous IR data. The mid-IR image 
shows a shell inside the well-known molecular borders of the nebula. The morphologies at 24 and 70 fim are quite different, and its 
color ratio is unusually warm. The far-IR image resembles the one at 8 /jm that enhances the structure of the molecular cloud and the 
Pillars of creation. We use this set of data to analyze the dust energetics and properties within this template for Galactic SFR. 
Methods. We measure IR SEDs across the entire nebula, both within the shell and the PDRs. We use the DUSTEM model to fit these 
SEDs and constrain dust temperature, dust size distribution, and interstellar radiation field (ISRF) intensity relative to that provided 
by the stai' cluster NGC661 1. 

Results. Within the PDRs, the dust temperature, the dust size distribution, and the ISRF intensity are in agreement with expectations. 
Within the shell, the dust is hotter (~ 70 K) and an ISRF larger than that provided by NGC661 1 is required. We quantify two solutions 
to this problem. (1) The size distribution of the dust in the shell is not that of interstellar dust. (2) The dust emission arises from a hot 
(~ 10* K) plasma where both UV and collisions with electrons contribute to the heating. 

Conclusions. We suggest two interpretations for the Ml 6s inner shell. (1) The shell matter is supplied by photo-evaporative flows 
arising from dense gas exposed to ionized radiation. The flows renew the shell matter as it is pushed out by the pressure from stellar 
winds. Within this scenario, we conclude that massive star fonning regions such as M16 have a major impact on the carbon dust 
size distribution. The grinding of the carbon dust could result from shattering in grain-grain collisions within shocks driven by the 
dynamical interaction between the stellar winds and the shell. (2) We also consider a more speculative scenario where the shell would 
be a supernova remnant. We would be witnessing a specific time in the evolution of the remnant where the plasma pressure and 
temperature would be such that the remnant cools through dust emission. 

Key words. 

1. Introduction 12004 )) and MIPS 24 //m ( iRieke et alj |2004|) . show a shell- 

^ ' like emission at 15 and 24 jum that fills the nebula cavity 

k'm ^""^^^ ^^^^^ "^"^^^ " - ''P''' llFlagevetalJl2009i) . as dehneated by the shorter IR wave- 

p><; LHdlenbrand et al. 1993) massive star formmg region made a sky ^^^^^^ the extent of the H, emission. The shell stands out 

H ■ icon by the publication of spectacular Hubble Space Telescope j„ j^e ISO 15 /urn and MIPS 24 jjm images, while the Nebula 

. (HST) images of the ionized gas emission ([ Hester et alj|1996l ). pjjj^.^^ the outer rim of the nebula are the strongest emis- 

As one of the nearest star forming region and one of the most gion features at 8 mn . Based on some spectroscopic evidence 

observed across the electromagnetic spectrum, the Eagle Nebula dUrguhart et al.ll2^ . we know that the mid-IR shell emission 

is a reference source. The nebula cavity is carved into the molec- ^^^^^ f^^^ ^^^t with only a minor contribution from ionized gas 

ular cloud by a cluster of 22 ionizing stars earlier than B3 jjnes to the broadband emission. 
(Dufton et al. 2006b') and with an estimated age of 1-3x10 yrs 

(Hillenbrand et al. 1993; Dufton etal. 2006b; iMartavan et alj M16 is not alone in this respect. There are other large, 

[20081) . partially symmetrical and rich HII regions (in terms of 

The mid-IR images of M16 either from the Infrared Space their OB stellar content) that displ ay a similar rnid-IR color 

Observ atory Camera (ISOCAM [C esarsky et al. 1996as!) at 8 and stratification: t he Rosette Nebula (iKraemer et alJ l2003h . the 



15 jum dPilbratt et al.lU998l; lOmont et al.li2003i) or base d on the Trifid Nebula ('Le flochet all 1 19991: iRho et alJ I2006'). and M17 



combined Spitzer observations using IRAC 8 pm (.Fazio et alJ (Povich et al. 2007^). Furthermore, the multi-wavelength obser- 

vations of the HII regions in the Galactic Plane, usin g the Spitzer 

Send offprint requests to: N. Flagey GLIMPSE and MIPSGAL Legacy surveys dChurchwell et al.l 
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l2009HCarev et alJl2009l) display overall a wide variety of com- 
plex morphologies, and show many "bubble"-like objects with 
a similar color stratification as M16 (Watson et al. 2008, 2009), 
although they are smaller and driven by one or a few OB stars. 

What are these Spitzer images of massive star form- 
ing regions teaching us about dust and the interaction of 
the stars with their environment? The IRAC and the MIPS 
24/im camera are imaging the emission from PAHs and Very 
Small Grains (VSGs). A first key to the interpretation of 
Spitzer images is the change in abundance and excitation 
of these small dust particles from molecular to ionized gas. 
Observations of nearby molecular clouds illuminated by O 
stars, where observations separate the H II photo-ionized 
gas layer from the neutral Photo-Dissociation Region (PDR) 
show that the PAH bands, which are a characteristic of PDR 
mid-IR emission spectra, are strikingly absent from that of 
the H II layer (e.g the Orion Bar and the M17SW interface, 
iGiardet al. 1994: Cesarskv et al. 1996b: Povich et al. 2007) . 
PAHs are quickly destroyed when matters flows across the 
ionization front. Several destruction mechanisms have been 
proposed: chemisputtering by protons and photo-thermo 
dissociation and/or Coulomb explosion associated with ab- 
sorption of high energy photons. Much less is known about 
the evolution of VSGs. The mid-IR shells may reflect dust 
processing by hard photons and shocks that impact the frac- 
tion of the dust mass in VSGs, but this possibiUty has yet to 
be constrained by modeUng of the dust emission. 

The evolutionary stage of the massive forming regions 
is a second key to the interpretation of the Spitzer im- 
ages. The mid-IR shells do not fit the classical view of the 
evolution of HII regions where the matter is swept away 
by the simultaneous effect of the ionization, stellar winds, 
and radiation pressure from their central OB stars (e.g. 
Tenorio-Tagle et al. 1982; Beflrametti et al. 1982; Rozyczka 
1985). In this scenario, the HII regions are "hollow". One in- 
teresting possibiUty is that gas photo-evaporating from dense 
condensations exposed to ionized radiation, creates a gas 
mass input within the cavity suflicient to balance the outward 
flow of matter. Are the shells reflecting such a mass input? To 
show that this is a plausible interpretation, one must quan- 
tify the mass input, as weU as the dust properties and excita- 
tion conditions, required to match the shells brightness and 
its distinct mid-IR colors. 

So far most of the studies on the mid-IR properties 
of these HII regions and smaller bubbles have been phe- 
nomenological and looking into the spatial distribution of the 
different emission components and not their physics. A small 
bubble where a more quantitative analyis has been carried 
out is G28.82-0.23 (Everett & ChurchweU. 2010). G28.82- 
0.23 (aka N49) is nearly sphericaUy symmetric, excited by a 
single 05V star, which has a thick 8 yum sheU surrounding at 
24 i^m a diffuse bubbl e (see e.g. Watson et al. 2008, Fig. 7). 
[Everett & Churchwell (2010) proposed a model where the 
mid-IR emission of G28.82-0.23 arises from dust entrained 
by the stellar wind. This interpretation involves a hot (> 10^ 
K), high pressure plama (p/k ~ 10^ K.cm^) where dust life- 
time is shorter than the expansion timescale. It requires that 
dust is constantly replenished by photo-evaporation of high 
density (10^ cm dusty gas cloudlets that have been overrun 
by the expanding nebula. CoUisional excitation by hot elec- 
trons contribute significantly to the heating of dust. Infrared 
dust emission is the dominant cooling channel of the dusty 
wind, which reduces the energy available for wind-driven ex- 
pansion. It seems to us that this specific model does not of- 



fer a general framework to interpret observations of larger 
HII regions, where one observes a similar 8 and 24 //m color 
stratification. 

The motivation of this paper is to study the nature of mid-IR 
shells in massive star forming regions using the Eagle Nebula 
as a template source. The detailed data available on this nearby 
nebula allow us to perform a quantitative modeling of the dust 
heating by UV radiation and, possibly, by collisions in a hot 
plasma. We quantify the dust emission in terms of dust physics, 
before discussing possible interpretations within an evolution- 
ary scenario of the Eagle Nebula as a massive star forming re- 
gion. In section 12] we present the Spitzer imaging observations 
of the Eagle Nebula from the MIPSGAL Galactic plane survey. 
Section[3]describes the morphology of M16 based on IR photo- 
metric and spectroscopic observations. We measure the spectral 
energy distribution (SEDs) across the entire nebula combining 
data from the ISO, MSX and Spitzer space missions. In sec- 
tions |4] and HJ we present exhaustive modeling of the dust 
properties. We first model the dust SEDs with UV heating 
only, and this sets constraints on the radiation field intensity 
and dust size distribution. Then we consider the possibiUty 
that the shell emission arises from a hot plasma where dust 
would be heated by colUsions with electrons. The reader not 
interested in the details of the modeUng can skip sections S] 
andHJ In section |6l we propose two scenarios of the present 
evolutionary state of the Eagle Nebula, which could account 
for the mid-IR shell and fit within present observational con- 
straints. The paper results are summarized in section Q 

2. Observations 

The Eagle Nebula has recently been observed by the Spitzer 
Space Te lescop e as p art of the GLIMPSE (program #00146, 
Beniaminltal] llool and MIPSGAL (program #205976, 
iCarev et al ■1 120091)" inner Galaxy surveys. The GLIM PSE survey 
has made use of the Infrared Array Camera (IRAC, iFazio et aU 
2004), while MIPSGAL has been reaUzed with the Multib and 
Imaging Photometer for Spitzer (MIPS, Rieke et al. 20 04|). In 
both cases we have used their enhanced products (S quires et al.l 
120051) . The MIPSGAL 24 /im data has been complemented with 
archival observations (Spitzer program #20726) and reprocessed 
using the standard Spitzer Post-Basic Calibrated Data toolfl A 
three-color image combining IRAC and MIPS data is shown on 
figure [1] 

Most of the data processing per formed on the MIPSGAL 
24 fjm o bservations is described in iMizuno et al ] (120081) and 
I Carey et al. (2009). At 70 pm, Spitzer detectors are Ge:Ga pho- 
toconductors. When observing bright, structured emission, like 
the one in the Eagle Nebula, such detectors show significant 
variations in responsivity, which manifest themselves as visible 
stripes in the final images, and result in photometric errors of 
several tens of percent. This effect has required an offline re- 
processing of the data, with tools specifically designed to, at 
the same time, reconstruct the history-dependent responsivity 
variations of the detectors and mitigate the associated stripes. 
The photometric uncertainty of extended emission is lowered 
from about 50% on the brightest features down to about 15% 
on the enhanced MIPS 70 pm data. The specific pipeline devel- 
oped for the MIPSGAL 70 pm observations will be detailed in 
Paladini et al. (in prep.). 

We complete the Spitzer observations of M16 with previous 
IR survey from MSX and observations from ISO, both photo- 



http://ssc.spitzer.caltech.edu/postbcd/ 
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Fig. 1. Composite Spitzer color image combining the IRAC 5.8 jum (blue) bands with MIPS 24 (green) and 70 fim (red). The FOV is 
~ 30', N is up and E is left. The two black boxes outline the Pillars of Creation, which raise from the bottom to the center, pointing 
slightly to the West and the Spire, on the East, almost pointing straight toward the West. The position and spectral type of the most 
massive stars of NGC661 1 is overplot: O stars are in red, B stars are in white. 



metric and spectroscopic . The ISOCAM/CVF spectra have al- 
ready been presented by lUrquhart et al.l (l2003l) . A slice of the 
ISOCAM/CVF spectroscopic cube is shown on Fig. |2(d)[ 



3. Observational results 

We use the many IR observations available to create a portrait of 
the nebula from NIR to FIR wavelengths. We then perform aper- 
ture measurements on both the broad band images and spectro- 
scopic observations in order to get characteristic spectral energy 
distributions (SEDs) and spectra of the Eagle Nebula. We focus 
our comments on the two main features of the nebula: the PDRs 
and the inner shell. 



3.1. Images 

The three-color image of Fig. [1] clearly highlights differences 
between intermediate wavelengths on the one side (MIPS24 in 
green) and the shorter and longer wavelengths on the other side 
(IRAC8 in red and MIPS70 in blue). The whole molecular cloud 
appears in purple while the inner shell is green. 

- At wavelengths shorter than ~ 10 yum, IRAC, MSX and ISO 
observations show the molecular cloud surface heated by the 
cluster UV radiation. The Pillars of Creation, the Spire (see 
Fig.[T]to identify these structures) and less contrasted emis- 
sion extend towards the cluster from the North and the East. 
To the NW and the SE, the rim of an outer shell can be iden- 
tified. It corresponds to the edge of the Eagle Nebula as seen 
in Ha. 

- At intermediate wavelengths, between ~ 12 and 24 //m, 
MSX, ISO and MIPS observations exhibit a significantly dis- 
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Fig. 2. ISOCAM/CVF mean spectra obsei-ved (a) on Pilbratt's 
blob, (b) at the tip of the main Pillar of Creation and (c) within 
the Pillar of Creations. Dotted lines are ON spectra, dashed lines 
are OFF spectra, thick solid lines are ON-OFF spectra. OFF and 
ON positions are shown on the ISOCAM/CVF 3' by 3' field of 
view, here at the wavelength of 12 jum. North is up. East to the 
left. 



tinct morphology, with a shell filling the inside cavity in be- 
tween the Pillars and the edges of the molecular cloud seen 
at shorter wavelengths. The shell extends over ~ 12' in the 
NW-SE direction towards the pillars and further out to the 
SW where there is no emission either at shorter or longer 
wavelength. There are some bright features within the shell, 
some of which have already been identified ( e.g. Pilbratt's 
blob, to the East of the main Pillar of Creation Pilbra tt et alJ 
|l_998). The lack of far-infrared observations prevented pre- 
vious authors to conclude anything specific on the nature of 
this shell. 

- At longer wavelength, the MIPS 70 fim observations are very 
similar to those at shorter wavelengths and mainly show the 
molecular cloud surface. The diffuse emission within the in- 
side cavity is visible but not as bright as at intermediate 
wavelengths. The lower angular resolution of these obser- 
vations does not allow us to make more detailed comments 
at this point. 

The IR morphology of the Eagle Nebula i s com mon among 
other star forming regions. Church well et al. (■2006^ have listed 
many such "bubbles" across the entire GLIMPSE Galactic plane 
survey with IRAC. Combining GLIMPSE and MIPSGAL 24 ^^m 
surveys reveals an inner shell for most of these regionsQ. 

3.2. SEDs measurements 

We perform ON-OFF aperture measurements to get both spec- 
troscopic and photometric SEDs. As shown on figure [2(d)] there 
is a band of unavailable pixels on the ISOCAM/CVF observa- 
tions. This band goes exactly through interesting and contrasted 

^ http://www.spitzer.cahech.edu/Media/releases/ssc2008- 
ll/ssc2008-lla.shtml 



features Uke the tip of the main Pillar and Pilbratt's blob. Rather 
than linearly interpolate the missing pixe ls like it has been done 
previously on ISOCAM/CVF data (e.g. lUrqu hart et al. 2003i 
we use these data as is. We present and interpret spectroscopic 
and photometric measurements separately. 

3.2.1. Spectroscopic measurements 

We compute average spectra on multiple positions within the 
Pillars of Creation area covered by the ISOCAM/CVF data. 
We use square boxes of 4x4 pixels (24x24"on ISOCAM/CVF 
6"pixel field of view) to estimate the mean brightness of sev- 
eral features. We use this method for both "ON" and "OFF" 
positions. We combine three different OFF positions to build a 
unique OFF spectrum. The resulting ON-OFF spectra are shown 
on figure |2] for two positions within the main Pillar of Creation 
and one on Pilbratt's blob. These three positions, marked on fig- 
ure [2(^]_corres£ondrespectively to spectra D, B and A of figure 
2 from lUrquhart et all (12003). One of our OFF positions is close 
to their spectrum C. As a consequence, our results are similar to 
theirs: 



- The spectra of the Pillars of Creation (see Fig. |2(b)| and |2(c)[ ) 
exhibit the characteristics of PDRs spectra with strong PAH 
features and gas lines. They also present the Si absorption 
feature around 10 fim. There are some variations between the 
two positions, mainly regarding PAHs features and gas lines 
strength, which traces variations in the excitation conditions 
between these two positions within the column of gas and 
dust. 

- The spectrum of Pilbratt's blob (see Fig. |2(a)) exhibits a 
strong continuum with very weak gas lines and PAHs bands. 
We thus assume, as a first approximation, that the MIPS 24 
yum shell is dust continuum dominated. 

- The OFF position has a spectrum with a weaker continuum 
than the blob but stronger than the Pillars. It also has much 
weaker lines and features than within the gaseous and dusty 
columns. 

3.2.2. Photometric measurements 

We combine IR observations of the Eagle Nebula from three 
different observatory: MSX, ISO and Spitzer Therefore, we 
first lower the spatial resolution of each observations to that of 
MSX data (20")- Then, as we did with the spectroscopic mea- 
surements, we pick up several interesting and contrasted fea- 
tures within the nebula. We name them as follows. The "PDR" 
group of features contains the tip of the main Pillar of Creation 
("Pillar", also known as Column I, with an embedded source at 
its tip, see Fig. |3(a)| i, the tip of the Spire ("Spire", also known 
as Column IV, with an embedded source at its tip, see Fig. |4(a)[ ) 
and a PDR within the main Pillar of Creation ("Shoulder", see 
Fig. |5(a)[ ). The "Shell" group of features contains Pilbratt's blob 
("Blob", see Fig. |6(a)| i, the contrasted border of the main shell 
("Shell border", see Fig. |7(a)) , a diffuse shell that extends to- 
wards the opposite direction ("Reverse shell", see Fig. |8(a)| i, a 
bright filament on the North- West side of the nebula("Filament", 
see Fig. |9(a)| i and some more diffuse emission on the South- West 
side of the nebula ("Diffuse", see Fig. |10(a)] i. For each structure, 
the main difficulty of the measurement is to properly estimate 
the background emission behind each of them. This is particu- 
larly true for the MIPS 70 fj.m images. 

We illustrate our method on the example of Pilbratt's blob 
but it is mainly valid for the whole set of structures. We first 
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Fig. 3. (a) Three color image as in figure [T] with the region Fig. 5. Same as figure|6]for the position of the "Shoulder", 
along which the profiles are measured for the main Pillar of 
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(b) select a rectangular area that encompasses the blob, as shown on 

Fig. |6(a)| We choose the orientation of the selected area in such 
Fig. 4. Same as figure|6]for the position of the "Spire". a way we avoid to select other neighboring contrasted features 

(e.g. the Pillars of Creation). We then compute the mean profile 
of the blob and its surrounding by averaging all the pixels along 
the short axis. The resulting normalized profiles for Pilbratt's 
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(b) 



Fig. 7. Same as figure |6] for the position of the "Shell Border". Fig. 9. Same as figure|6]for the position of the "Filament" 
The darker sections of the profiles show the top and bottom of 
the "jump" used to measure the fluxes at each wavelength. 
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(b) 



Fig. 10. Same as figure|6]for the position of the "Diffuse" 



Fig. 8. Same as figure |6]for the position of the "Reverse Shell". 

We then measure the mid to far-IR SED of each structure. 
We adapt the method as a function of the profile shape. For the 
structures that present a peak of emission at every wavelength 
blob are shown on Fig. |6(b)| for several wavelengths. The profiles (e.g. Pilbratt's blob. Spire), we estimate the background through 
for the other features are shown on Fig. |7(b)| to |10(b)| a spline interpolation of the profile on both sides of the peak 
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(see Fig. |5(b)[ l. The flux of the structure is thus given by the inte- 
gration of the background subtracted profile over the size of the 
structure. The actual size over which we integrates the flux may 
slightly vary from one channel to another. The uncertainty on 
each measurement is given by the range of background values 
as estimated by the spline interpolation. For the other structures, 
where the profiles exhibit a "jump" (case of the shell border, see 
Fig. |7(b)| i, we estimate the height of the "jump" at each wave- 
length by measuring the difference of the surface brightness be- 
tween the top and bottom of the "jump". The uncertainty on each 
measurement is given by the standard deviation of the surface 
brightness at the top and the bottom of the "jump". 

While the measurements are usually straightforward on the 
MIPS 24 yum profiles, they are significantly more uncertain on 
the MIPS 70 fim profiles, especially for less contrasted struc- 
tures like the "Filament" or the "Diffuse" emission. In those two 
last cases, we are not sure about the exact spatial extent of the 
structure at 70 fim and the range over which to estimate the back- 
ground (see Fig. |10(b)| i. This generally also applies to the IRAC 
8 and 6 fim measurements, but to a lesser extent. In particular, 
for the "Filament" structure, the discrepancy in the profile's peak 
position between MIPS 24 /urn and MIPS 70 //m or IRAC 8 yum 
is significant enough so we do not consider them as probing the 
same physical conditions (see Fig. |9(b)| and |10(b)] i. Since there 
is no other obvious feature at the position of the MIPS 24 yum 
peak, we will thus use the MIPS70 i-im measurement as an upper 
limit. Additionally, the uncertainty on the MIPS 70 fim flux of 
the "Diffuse" is significantly higher. The resulting photometric 
SEDs, normalized to MIPS 24 fim, are presented on figure [TT] 
Again, the differences between the structures within the shell and 
those within the PDRs are clear. 

- The PDRs of Ml 6, both at the tip of the Spire and within the 
Pillars of Creation, are characterized by an almost flat SED 
from near to mid infrared and a continuous increase mid to 
far infrared wavelengths. The SEDs of the position with an 
embedded source ("Pillar" and "Spire") do not appear to be 
different from that of the "Shoulder" at near infrared wave- 
lengths. At longer wavelengths, the SED of the "Shoulder" 
increases slightly less than those of the "Pillar" and the 
"Spire", which encompass embedded source. The ratio be- 
tween MIPS24 and MIPS70 is about 0.1 for the "Shoulder" 
and about 0.3 at the tip of the main Pillar of Creation and the 
Spire. 

- The inside shell, at Pilbratt's blob position and on blight con- 
trasted structures, is characterized by a significantly steeper 
increase of the intensity from near to mid infrared and a flat 
or decreasing SED from mid to far infrared. On Pilbratt's 
blob, the Shell border and the Reverse shell, the MIPS24 to 
MIPS 70 ratio is about 4.5, 2.3 and 0.95 respectively. 

- The Filament and the Diffuse SEDs appear in between these 
two sets of SEDs. Both their MIPS24 to MIPS 70 ratio is 
lower than inside the shell and their near to mid infrared SED 
is steeper than within PDRs but the uncertainties are signif- 
icantly larger As a consequence, in the following sections, 
we do not discuss these last two positions. 

The measurements of the near-IR to far-IR SEDs con- 
firm what spectroscopic observations were suggesting: the dust 
within the inner shell is significantly different from that within 
PDRs. The addition of the MIPS 70 fim and its comparison to 
MIPS 24 fj.m provide us with constraints on the position of the 
dust emission peak in the FIR. We explore in the next section 
whether the difference arises from external excitation or intrin- 
sic properties. 
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Fig. 11. Comparison of the structures SED. Solid lines: struc- 
tures within the shell. Dotted lines: structures within the PDRs. 
Red dash: Filament. Green dash: Diffuse 

4. UV heating of the dust 

In this section , we model the du s t emis sion within M16 using the 
dust model of'Compieg ne et al] (1201 lb . In this model, the dust is 
heated by the incident flux of UV photons only. We first show 
that the MIPS 24 jjm to MIPS 70 yum ratio may be direcdy re- 
lated to the intensity of the interstellar radiation field (ISRF) in 
the shell. We then use the dust model to determine the best set 
of parameters that describes the complete observed SEDs over 
the entire nebula. In this section, we limit ourselves to the fol- 
lowing parameters : the intensity of the incident radiation field 
and the dust size distribution, in terms of abundance of the dust 
components. 

4.1. Method 



The dust model of 'Compiegne et al.' ("201 1") is an updated ver- 
sion of the original Desert et al. (1990) model. In their model, 
iCompiegne et al.l (1201 1 ) use four dust components: (1) poly- 
cyclic aromatic hydrocarbons (PAH), (2) stochastically heated 
very small grains of amorphous carbon (VSG or SamC), (3) 
large amorphous carbon grains (LamC) and (4) amorphous sili- 
cates (aSil). We combine LamC and aSil grains into a unique big 
grains (BG) component using these grains relative abundances 
found in the diffuse hi gh galactic latitude (DHGL) medium 
jCompiegne et al.ll201 Ih . We assume a fix dust-to-gas mass ra- 
tio of 1%. We then use the dust model to compute the emission 
spectra of the three dust components (PAHs, VSGs and BGs) 
illuminated by the incident radiation field from the star cluster 
NGC6611. 

We use fli e STA RBU RST99 onhne modeP I descr ibed in 
Leit herer et al ] (Il999h and IVazquez & Leithere J(l2005l) to de- 
fine the spectral shape of the radiation field from the illumi- 
nating star cluster NGC6611. We use the following parame- 
ters: 2 millions years old cluster, Salpeter initial mass function 
(dn/dM oc M'^-^^), stellar masses from 1 Mq to 100 Mq. The 
modeled radiation field coiTesponds to 1 .6 x 10'' L©. We normal- 
ize it so that it is in agree ment with the total flux of the most 
massive stars of the cluster. IPufton et an (l2006b) have presented 
an analysis of VLT- FLAMES spectroscopy for NGC661 1. Their 
online catalogue (Dufton et al. 2006a) lists stars classified as ear- 
lier than B9. The 42 members of NGC661 1 have a combined to- 
tal luminosity of 3.4 x 10^ Lq, which is a factor 480 smaller than 



http://www.stsci.edu/science/starburst99/ 
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the Starburst99 model output spectrum. We apply that correction 
factor to the model spectrum of the ISRF. In Habing units - in- 
tegrated intensity of the solar neighborhood from 912 to 2000 A 
or 1.6 X 10"^ erg.s^'.cm"^ - the cluster radiation field intensity 
is xo - 4800 at a distance of 3 parsecs (see section l4~4l for a dis- 
cussion on the spatial variations of the IRSF). In the following, 
we use this value as a reference for the dust model. 

For the features within the shell ("Blob", "Shell border" and 
"Reverse shell"), the use of a non-attenuated radiation field is 
acceptable since the UV optical depth is low. For the features 
within the PDRs ("Pillar", "Spire" and "Shoulder"), we have to 
take into account the extinction of the ISRF by the ionized layer 
of gas and the FDR layer itself. We model this in a simple way 
by removing the Lyman continuum photons and with a far-UV 
extinction of 1 magnitude. Such an extinction accounts for the 
fact that the emission from PDRs comes from a range of depths 
into UV-dark clouds with a weighting proportional to the UV 
field. A more detailed study of the PDRs is beyond the scope of 
this paper. 

4.2. MIPS 24 fim to MIPS 70 urn ratio as a tracer ofx 

We first use the dust model of ICompiegne et all (1201 Ih to com- 
pute the MIPS 24 jum to MIPS 70 //m ratio of the dust emis- 
sion for different dust size distributions to show how it is re- 
lated to X- Within this wavelength range, the PAHs contiibu- 
tion to the emission is weak relative to that of VSGs and BGs. 
Therefore, we present the MIPS 24 jum to MIPS 70 yum ratio 
as a function of x for three size distributions: VSGs only, BGs 
only and a mixture of VSGs and BGs that matches their relative 
abundance in the diffuse high Galactic latitude medium (DHGL, 
ICompiegne et al.l [201 1,) . Therefore, we take into account any 
dust evolutionary process that would destroy a specific grain size 
component. Figure [12] shows the results along with the MIPS 24 
/Ym to MIPS 70 yum ratio measured for the Eagle Nebula struc- 
tures, both within the shell and the PDRs. The differences be- 
tween the set of curves for the PDRs and that for the shell are not 
significant. We first make no distinction while presenting them. 
Then we discuss the results for the PDRs and Shell structures 
independently. 

For a given xlx^'iy VSGs always have a higher 
MIPS24/MIPS70 as they are hotter than BGs. However, 
forxIXQ ^ 1-0' MIPS24/MIPS70 is almost independent, within 
a factor of a few, from the dust size distribution. These values 
of X correspond to the large values of the MIPS24/MIPS70 
(> 1). For xlXo ^ 1-0, MIPS24/MIPS70 is significantly 
more dependent on the grain size distribution with difference 
up to almost 2 orders of magnitude. Likewise, for a given 
MIPS24/MIPS70, the required xlxo is always higher for BGs 
than VSGs. The difference is as small as a factor of a few for 
high values of MIPS24/MIPS70 and as high as almost 2 orders 
of magnitude for low values of MIPS24/MIPS70. Therefore, 
given MIPS24/MIPS70, the constraint on the intensity of the 
IRSF is stronger for higher values of x and requires a better 
knowledge of the dust size distribution (e.g. as provided by 
other IR observations, see next subsection) at low values of x- 
On the contrary, constraining the dust size distribution requires 
an a priori onx and can better be done at low values of;^'. 

According to the model, the PDR stiTictures ("Pillar", 
"Spire" and "Shoulder") require an ISRF intensity at most a fac- 
tor 2 lower than the reference, and no lower limit can be esti- 
mated because we have no constraint on the dust size distribu- 
tion. However, if we assume it does not significantly depart from 



that of the DHGL medium, the MIPS24 to MIPS70 ratio within 
PDRs are best interpreted xlXiSRFfi - 0.1. The inner shell 
structures ("Blob", "Shell border" and "Reverse shell") are more 
on the high end of the ISRF intensity. The "Shell border" and the 
"Blob" are best inteipreted with ^xlXiSRFf.) of at least a few and 
up to 16, whether the dust size distribution is dominated by BGs 
or VSGs. The difference between the ISRF intensity that illu- 
minates these two structures and the PDRs is thus at least an 
order of magnitude. The "Reverse shell"' position however is 
not strongly constrained and overlaps those of the PDRs struc- 
tures. If at this position the dust size distribution is dominated 
by VSGs, ih&nxIxisRFfi - 0.1 ^^\\& xIxisrf,() ~ 1 if the BGs 
contribute the most to the dust size distribution. The full range 
of required ISRF intensities for each structure is given in Table 

m 

Table 1. Lower and upper limits of xlxo for the whole set of 
structures as deduced from their MIPS 24 yum to MIPS 70 jjm 
ratio. 



Shell structure 


XlXo 


PDR structure 


XlXo 


Reverse shell 


0.13-1.3 


Pillar 


< 0.3 


Blob 


5.6-16 


Shoulder 


<0.6 


Shell border 


1.4-5.4 


Spire 


< 0.4 



Indirectly the MIPS 24 to MIPS 70 //m ratio also pro- 
vides us with a measurement of the equilibrium dust tempera- 
ture Teq of the largest dust particles. In figure [13] we plot the 
BGs equilibrium temperature, provided by the dust model, as a 
function oix^ for both the PDR and the Shell s tructures, and for 
both types of large grains used in the model of Compiegne et al.l 
(2011): LamC and aSil. For a given radiation field intensity 
X/XisRF.o, we plot the upper and lower limits for the equilib- 
rium temperatures of each grain type. The difference between 
both type of BG components is not really significant. In figure 
[T3l we hatch the range of equilibiium temperatures for the val- 
ues of ;i'/;^'o given by Fig. [12] 0.13 <x/Xo < 16 for the Shell and 
X/XQ < 0.6 for the PDR structures. While the smallest LamC 
grains in the PDR structures may reach equilibrium temperature 
as high as 71 K, those are in limited number Likewise, only the 
largest grains in the Shell structures may reach equilibrium tem- 
perature as low as 24 K. The majority of the grains, as traced by 
the most abundant size bin of each BG component (also plotted 
in figure [T3Tl. span a range of equilibrium temperatures that does 
not overlap significantly between the Shell and the PDR struc- 
tures. For the PDRs structures, equilibrium temperatures for the 
most abundant size 20 < Tgg < 50K while for the inner shell 
structures 35K < T^q < IQQK. Therefore, equilibrium temper- 
atures above 50 K can only be efficiently reached by BGs in 
the inner shell while equilibrium temperatures below 50 K are 
mostly found in the PDRs. The dust in the inner shell is thus sig- 
nificantly hotter than that in the PDRs. Ilndebetouw et alj (|2007|) 
have used IRAS 60 yum to IRAS 100 yum ratio to build a low spa- 
tial resolution (4.3') color temperature map of the dust in M16. 
Their values range from 32 K in the molecular cloud to 40 K in- 
side the nebula. We build the same map (not shown here) using 
IRAS 25 yum to IRAS 60 yum ratio (to better match the MIPS24 
to MIPS 70 yum diagnostic) and find color temperature ranging 
from 45 K to 65 K, more in agreement with our measurements 
of the BGs equilibrium temperature in the shell. The remaining 
difference may come from the lower spatial resolution that aver- 
ages "hot" features with "cold" features within the beam. 
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Fig. 12. MIPS 24 fim to MIPS 70 fim ratio as a function of the ISRF intensity, as predicted by the model oflCompi egne et al] ( 1201 lb . 
Several dust size distribution are used: (dotted line) BGs only, (dashed-line) VSGs only and (solid line) mix of BGs and VSGs. The 
MIPS24-to-MIPS70 ratio for several structures within M16 is indicated. The ISRF spectral shape is that mention in the text with (a) 
no extinction, (b) A(FUV) = 1 mag and the Lyman continuum photons removed. 
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Fig. 13. BGs equilibrium temperature as a function of the ISRF intensity. The hatched area corresponds to the range of equilibrium 
temperatures span by the entire BGs size distribution. The solid lines represent the equilibrium temperature for the most abundant 
size bin. The hatched area and the solid lin e are only plotted for the values of x/xo that are given by figure [12] Black is for LamC 
grains, red is for aSil grains as described in ICompiegne et al.l (|20I ll) . The ISRF spectral shape is that mention in the text with (a) no 
extinction, (b) A(FUV) = 1 mag and the Lyman continuum photons removed. 



Table 2. Best-fit parameters for SEDs of the Eagle Nebula. The ISRF intensity, the dust size distribution, in terms of relative 
mass ratio a bundances, and the total dust column density are given. The parameters for the diffuse high Galactic latitude (DHGL) 
reference of lCompiegne et aD(l20I II) are also given. The dust-to-gas mass ratio is fixed at 0.01 therefore a dust mass column density 
of \ .1 ng .cnT^ corresponds to 10^" H.cm 



Position 


X/Xo 




Yvsa(M/MM) 


Ybg{M/Mh) 


cta,,,, (yug.cm ^) 


DHGL 




7.8 X 10-* 


1.65 X IQ-* 


9.25 X 10-' 


1.7 


Pillar 


0.19 ±0.04 


(2.64 ± 0.57) X 10-" 


(2.45 ± 0.90) X 10-" 


(9.49 ± 1.82) X 10-" 


380 


Shoulder 


0.43 ± 0.08 


(2.51 ± 0.45) X 10-'' 


(1.12 ± 095) X 10-" 


(9.64 ±2.15) X 10-3 


33 


Spire 


0.12 + 0.05 


(2.96 ± 1.27) X 10-'' 


(5.09 ± 2.89) X 10-" 


(9.20 ± 3.62) X 10-3 


870 


Shell border 


4.36 ± 1.36 


(4.85 + 1.12) X 10-" 


(3.69 ± 2.71) X 10-" 


(9.63 ± 2.77) X 10-' 


02 


Blob 


9.69 + 2.33 





(5.98 ± 3.07) X 10-" 


(9.40 ± 1.82) X 10-3 


2.9 


Reverse shell 


1.15 ±ai3 





(1.99 ± 031) X 10-3 


(8.01 ± 023) X 10-3 


2.1 


Shell Border 


2* 


(6.68 + 4.47) X 10-' 


(1.05 ± 077) X 10-- 


(3.59 ± 2.29) X 10-" 


017 


a(i{VSG) = 5.5 nm 
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Fig. 14. Best-fit for (a) the "Pillar", (b) the "Shoulder", (c) the "Spire", (d) the "Shell border", (e) the "Blob", (f) the "Reverse shell". 
Solid line is the total model spectrum, dotted-lines are PAHs, VSGs, and BGs contributions. Diamonds are model broadband fluxes. 
Red crosses are measurements. 



4.3. F/ff/ngf of \he whole IR SED 

The additional measurement provided by MIPS 70 fim enables 
us to give some constraint on the ISRF intensity that is required 
to heat the dust up to the observed temperatures. Hereafter we 
use our dust model and the whole IR SED of each structure 
within M16 to better determine the variation of x and the dust 
size distribution at the same time. 

We set five parameters free: the intensity of the ISRF and 
the abundances of the three dust components (PAH, VSG and 
BG). The spectral shape of the ISRF is that described in section 
14.11 The other parameters describing the dust size distribution 
(e.g. the size range and d istribution shape) are those presented in 
ICompiegne et all (l201 ih . We use the MPFIT packag^ for IDL 
(iMarkwardtl l2009h to constrain the free parameters, given the 
SED. We use the default tolerance parameters and limit the four 
parameters to positive values. The best-fit spectra are shown on 
Fig. [T4]and the best-fit parameters are given in table |2] We also 
give an estimate of the dust column density for each feature as- 
suming a dust-to-gas mass ratio of 0.01. 

The positions within the PDRs are best fit with low values 
of x/xo (a fsw 10"') in agreement with the upper- limits from 
Table [T] a factor of a few less PAHs and a factor of a few more 
or less VSGs than in the DHGL medium. An increase/decrease 
of the small grains abundance by a factor of a few within PDRs 
of NGC2 023N and the Horsehead Nebula has already been ob- 
served by ICompiegne et al.l (l2008l) and within translucent sec - 
tions of the Taurus Molecular Complex bv lFlagev et al.l (l2009bl) . 
The low values of x/Xo required to fit the SED of the PDRs 
can partly be explained with shadow efifects within the neb- 
ula. Another parameter that we do not take into account in our 
simple model is the geometry of the features and the resulting 
limb brightening effect. Indeed, the dust column density for the 
"Pillar" and the "Spire" position is about a few 10 g.cm - 



http ://purl . com/net/mpfi t 



which corresponds to a gas column density of a few 10^^ cm"^ 
or a visual extinction of a few magnitudes, significantly larger 
than that required for attenuating the incident UV radiation field. 
The three "PDR" positions give very similar results, especially 
in terms of PAH abundance which varies by less than 10%. The 
VSG abundance is varying more significantly, up to a factor 5. 
The BG always dominates the dust size distribution with abun- 
dance very close to that of the DHGL. 

The positions within the "Shell" require larger values of^/xo 
(about a few) in agreement with values from Table [Tj a signifi- 
cant depletion of the PAHs and a significant increase of the VSG 
abundance, up to a factor 10, with respect to the PDRs values 
and at the expense of the BG component. The total dust col- 
umn density is about ~ 10"^ g/cnr, similar to DHGL values and 
which corresponds to a gas column density of about 10^° cm"^. 
As a consequence of the increased the VSG and BG emission 
spectra peak at very close wavelengths (see Fig. |14(d)]|14(e)| and 
|14(f)| i. We show in the previous section that MIPS24/MIPS70 is 
a good tracer of but not of the dust size distribution, especially 
at high values of x- Here, the addition of the other IR obser- 
vations provides better constraints on the VSGs to BGs relative 
abundance. For the position of the "Reverse shell", the initial 
best-fit (not shown here) underestimates the MIPS 70 //m mea- 
surement by almost an order of magnitude. As a consequence, 
the required is overestimated relative to that from Table [1] 
derived from the MIPS24 to MIPS70 ratio. We believe this poor 
fit at the longer wavelengths is due to the uneven number of mea- 
surement at short and long wavelengths, relative to the peak of 
the dust emission. From 6 to 24 fim, no less than seven measure- 
ments are available while only MIPS 70 fim is available at wave- 
lengths longer than the peak position. The fit process is thus bi- 
ased towards shorter wavelengths. In order to limit this effect, we 
repeat the fitting process of the "Reverse Shell" position with an 
increased weight on the MIPS 70 i-im measurement. Figure [l"4(Q] 
shows the result of that fit. The three positions within the shell 
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give results that are very similar to each other and very different 
from those of the PDRs positions: (1) an incident radiation field 
intensity a factor of a few larger than that provided by the star 
cluster NGC6611 and about an order of magnitude larger than 
that required for the PDR positions, (2) a significant depletion of 
the PAHs and (3) an increase of the VSGs abundance relative to 
BGs as compared to the PDRs positions. 

In order to explore furthermore the importance of a 
change in the dust size distribution, we redo the fit of the 
"Shell border" with a fixed intensity of the radiation field 
X/xo - 2 an d a free mean size of th e VSG component (ao). In 
the model of Compieane et al. (2011) for the DHGL medium, 
the VSGs size distribution is assumed to have a log-normal 
distribution (with the centre radius ao -2 nm and the width 
of the distribution cr = 0.35 nm). We keep the width of the 
log-normal distribution constant and set free the centre size 
ao between 0.6 and 20 nm. The other free parameters for that 
fit are the abundances of the dust components, as previously. 
The best-fit is plotted in Figure [15] and the parameters are 
given in Table |2j A significant increase of the mean size of 
the VSGs, by almost a factor 3, is required. There are almost 
no PAHs, as in the previous fits. The BGs are about a factor 
3 less abundant than in the previous fit and about a factor 30 
less abundant than in the DHGL. The abundance of VSGs 
is about 60 times higher than in the DHGL medium, though 
the uncertainty remains large (~ 75%). Therefore, the "Shell 
border" SED requires that most of the dust mass is concen- 
trated into the VSGs component. Despite those variations of 
the dust size distribution, the total dust column density re- 
mains very similar to that of the fit with a fixed mean size for 
VSGs (0.17 instead of 0.20 //g em 2). We also try the same fit 
with^/;t'o = 1 but find that the uncertainties on the parame- 
ters are then significantly higher (> 100%). 

We conclude that the MIR shell SED can either be ac- 
counted for a significant change in the dust size distribution 
or by an additional source of heating besides the star cluster 
radiation field. In the following, we first discuss two sources 
of UV heating that may account for the values otx/xo > 1 
required to fit the "Shell" SEDs. The first one is related to 
the spatial variations of x due to the exact positions of the 
OB stars in the sky. The second originates in the Lyman a 
photons emitted by the hydrogen and absorbed by the dust 
grains. We then consider, in the next section, another heating 
process originating from colUsions with the gas. 

4.4. Spatial variations of fhe incident radiation fieid 

Depending on the exact positions of the main OB stars of 
NGC6611 within the Eagle Nebula, the local incident radiation 
field intensity may vary and thus explain the required values of 
X/xo- For the"cold" PDRs features, it is easy to explain values 
of x/Xo < 1 as the stars are not all together on the plane of 
the sky, additionally to probable shadow effects already men- 
tioned. However, the required values ofx/xo > 1 for the "Shell" 
structures cannot be accounted for by the same interpretation. 
In figure [T] we indicate the position and the spectral type of the 
members of NGC6611, according to Dufton et al. (20061). We 
compute the variations of the ISRF intensity xo as a function of 
the position, taking into account the luminosity and position of 
each individual member of the cluster We assume that all the 
stars and the "Shell" structures are in the same plane of the sky. 
Therefore, the values of the local ISRF intensity we compute 
are thus upper-limits and the corrected values of x/XQ required 
for the best-fits are lower-limits. All these values are reporte in 




Fig. 15. Same as Figure [l4(d)| but with;^/;^'o = 2 and a free mean 
size of the VSG component. 




Fig. 16. Blow out of the MIPS 24 jjm image at the position of the 
Blob. The position and spectral type of O stars from NGC6611 
are also reported. The red dashed circle, centered on the 08.5V 
star has a 26 arcsec radius (about 0.25 pc at the distance of Ml 6). 



Table |3] The corrections factors are about a factor of a few at 
most. The required values of x/Xo for the Shell Border and the 
Blob are still at least a factor 2 to 3 higher than that provided by 
the star cluster 

The position of the members of NGC66 1 1 also reveals that 
Pilbratt's Blob is very close to an 08.5V star, as shown on 
Fig. [16] This suggest a possible local action of the winds from 
this star. The shock provided by the winds may account for a 
local enhancement of the density within the shell and possibly 
for dust processing. The same interpretation does not hold for 
the "Shell border" and the "Reverse shell" position which both 



Table 3. Correction factors on;^'o from the dispersion of the stars 
in the sky plane and corrected ;i'/;^'o required for the best-fits. 



Position 


Correction 


Corrected ^/^o 




factor 


(best fit) 


Shell border 


< 1.5 


> 2.9 


Blob 


<4.5 


> 2.1 


Reverse shell 


< 6.8 


> 0.2 
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are away from any OB star, as also shown in Fig. [16] We discuss 
collisional heating in section |5] 

4.5. Lyman alpha photons heating 

We show here that Lyman a photons are not a significant heat 
source for the shell. Every Lyman a photons emitted by an hy- 
drogen atom, after multiple absorption and reemission by other 
hydrogen atoms, either succeed to escape the medium or is ab- 
sorbed by a dust grain. The Lyman a contribution to the dust IR 
brightness is 5 LvQ. = / neXnH+xa2Xhviyadl - EMxa2Xhviya, 
where EM is the emission measure and 02 the hydrogen recom- 
bination coefficient to levels 2 and higher. The equation assumes 
that all recombinations from excited levels produce a Lyo- pho- 
ton that is absorbed by dust. 

We compute the EM from Bry observations of M16 ob- 
tained at the Canada-France-Hawaii Telescope (CFHT). These 
observations will be presented in a future paper. They do not 
show a counterpart of the "Blob", but there is an increase of the 
Bry emission associated with "Shell border" of EM = 3.5 x 
10^ pc.cm The Lya photons total flux that we estimate from 
these measurements is S lya - 0.048 erg.s^'.cm"^. In compar- 
ison, the 24 yum brightness of the "Shell border" is 230 MJy/sr 
which corresponds to a bolometric intensity of 0.37 erg.s"' .cm"^ 
that we measure on the best fit (see Fig. |7(b)| i between 1 and 
1000 ;um. The extra heating provided by the Lya photons is thus 
about a factor 8 too small. 



5. Collisional heating of dust 

In this section, we face the difliculty of explaining the shell in- 
frared colors with UV heating by considering the possibility that 
gas-grains collisions provide additional dust heating. We quan- 
tify the conditions that would be required to fit the shell SED 
with a combination of ra diative + colli sional heating of dust. 

We use the work of iDwekl (Il987l) to quantify the heat de- 
posited in the grain by collisions with electrons as a function of 
grain size and plasma temperature. Like in section |4] we use the 
DUSTEM mod el with a combination of silicates and amorphous 
carbon grains (ICompiegne etal.ll2011 l). Since the DUSTEM 
code does not include collisional excitation, we wrote a spe- 
cific module to compute the distribution of grain temperatures 
for stochastic heating by both photons and collisions. This code 
takes into account the Maxwellian distribution of the electrons 
kinetic energy. The results of our calculations are illustrated in 
Fig. [T7]for carbon grains. The Spitzer colors Iv(8//m)/Iv(24yum) 
and Iv(24yL/m)/Iv(70/im) are plotted versus grain size for radia- 
tive heating by the mean Eagle Nebula radiation field, and radia- 
tive +collisional heating for a range of electron densities n^. The 
temperature of the electrons Te is fixed to 10^ K. Our specific 
choice of is not critical, because the colors depend mainly on 
the plasma pressure, i.e. the product n^ x Tg. Collisional heating 
has a significant impact on the infrared colors for pressures p/k 
larger than a few 10^ K.cm"^. The figure shows that both colors 
may be fit for pressures pjk - 1.9 xn^Tg ~ 5 x\Q^ K.cm"-' and 
a characteristic grain size of ~ 10 nm. For this plasma pressure, 
collisions with electrons dominate the heating of small grains 
with radii < 10 nm, while radiation is the main heating source 
for larger grains. To illustrate the ability of the dust model to 
fit the shell SED, we use a dust size distribution that combines 
a log-normal size distribution for very small carbon grains plus 
a power-law size distribution for silicates. We keep the relative 
fractions of dust mass in carbon grains and silicates to their in- 



terstellar values: 1/3 and 2/3, respectively. In Fig. [18] we show 
a fit of the "Shell border" SED obtained for ng - 30 cm"^ and 

= 10^ K. For this fit, the characteristic radius (i.e. the mean 
value of the log-normal size distribution) of the carbon VSGs 
is 6.5 nm. This value is somewhat smaller than the value that 
may be infeiTed from Fig. [17] because the silicates contribute to 
about half of the 70 jum flux. The figure shows that for a given 
plasma temperature the characteristic grain size is tightly con- 
strained by the Iv(8//m)/Iv(24yum) ratio. It depends on the plasma 
temperature because this constraint is related to the stochastic 
heating of the smallest grains by collisions with electrons. The 
model also allows us to estimate the dust mass in the shell. The 
dust surface density is 2 x lO^-'Mopc"^. Scaling this value by 
the full extent of the shell (4 pc radius), we find a total dust mass 
of3xlO-2 Mo. 

The pressure inferred from the modeling of the collisional 
heating may be compared with independent constraints on the 
pressure within the Eagle nebula. This comparison raises diffi- 
culties with, but does not fully rule out, the collisional heating 
solution. The gas pressure inferred from Hubble observations 
optical line emission from the faint end of the photo-evaporation 
flows arising from Pillar I is pJ k ~ 10^ K.cm"-' (see Fig. 7b, 
abscissa in lHester et al. This value sets an upper limit 

on the ambient pressure around the flows, which is lower than 
the pressure required for the collisional heating solution. One 
possible way out of this problem is that Pillar I is not embed- 
ded in the shell. The shell pressure can also be estimated from 
Pilbratt's blob. The blob is close to an 08.5V star known to be 
associated with the ionizing cluster of the Eagle Nebula (see 
Fig. [T6] l. Its morphology and position on one side of the star 
suggests that it traces a bow shock cr eated by a supersonic m o- 
tion between the shell and the star (Ivan Buren et al.lll990l) . If 
this interpretation is right, it sets a constraint on the shell pres- 
sure. At the standoff distance dg, i.e. the distance between the 
star and the edge of the blob, there is a pressure equilibrium 
between the wind pressure and the ambient pressure plus the 
ram pressure associated with the star motion. Hence, the wind 
pressure at the standoff distance, p^ - M^. x VkHAtt x d^), is 
an upper limit on the ambient pressure. From the 24yum image, 
do - 0.2 pc. We use the empirical relation between wind mo- 
mentum and stellar luminosity (.Kudritzki & Puis 2000) : for an 
08.5V star M,. x V,, ~ 2 x 10"^ Mo-yr"' x 10^ km.s"'. Hence, 
we find pn lk - 2 x 10^ K.cm"^, a value more than one order 
of magnitude smaller than the pressure required for the colli- 
sional heating solution. Here, the plausible way out would be 
that Pilbratt's blob is not a bow-shock. 

6. The nature of the mid-IR shell 

In this final part of the paper, we discuss the results from our 
dust modeling in the context of the Eagle Nebula massive star 
forming region. We have shown in the previous sections that 
the dust SED of the MIR shell cannot be accounted for by 
standard models (i.e. interstellar dust heated by UV radia- 
tion). We find two possible explanations. (1) The fraction of 
the dust mass in stochastically heated VSGs is much larger 
in the shell than in the diffuse intertsellar medium. (2) There 
is an additional source of heating which could be collisional 
heating in a high pressure plasma. Here we present two sce- 
nari that can explain either or both of these requirements. In 
the first one the mid-IR shell is a windblown shell, where the 
dust is heated by UV photons and where large grains have 
been ground into stochastically heated small particles. In the 
second scenario we investigate a more speculative hypothesis 
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Fig. 17. Spitzer colors I^(8yum)/Iv(24//m) and Iv(24/^m)/Iv(70/im) 
for carbon grains versus grain size. The solid line give the colors 
for radiative heating for the Eagle Nebula ISRF. The other lines 
show the impact of coUisional heating for a range of plasma pres- 
sure and a fixed temperature Tg of 10^ K. A good fit of the shell 
SED is obtained for x T,, ~ 3 x 10^ K cm"^ (see Fig.fTSl) 



where the shell would be a supernova remnant that would be 
cooling through IR dust emission. 

6.1. A wind blown shell 

In this first scenario, matter outflowing from dense condensa- 
tions and exposed to ionizing radiation from the stellar cluster, 
in particular the Eagle pillars, supply the shell with a continuous 
inflow of gas and dust. The mechanical pressure from the stellar 
winds push this matter outward, but the shell persists provided 
that its outward expansion is compensated by continuing photo- 
evaporation. Since the shell is within the ionizing boundary of 
the nebula, the diffuse matter in the shell is fully ionized. The 
gas density and column density are too small to absorb all of the 
ionizing radiation. To quantify this scenario, we apply the em- 
pirical relation between wind momentum and stellar luminosity 
(Kudritzki & Puis 2000) to each of the O stars in the cluster For 
a shell inner radius of 3 pc, we find that the winds pressure is 
Pmnds/k = 5 X 10^ K.cm"-'. This value is a few times larger than 
the radiation pressure estimated from the shell infrared bright- 
ness Pradlk ~ Bjr/c ~ 10^ K.cm""', wherc Bjr is the mean bolo- 
metric IR brightness ~ 0.4 erg.cm^^.s ' and c the speed of light. 
The shell matter moves outward, because the wind pressure is 
higher than the average pressure in the interstellar medium. The 
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Fig. 18. Fit the spectral energy distribution measured on the 
Eagle shell with radiative plus collisional heating. The ISRF is 
that determined in section|4]with;^/;t'o = 1. The electron density 
is 30 cm"^ and the plasma temperature 10^ K. 



expansion velocity is commensurate with the sound speed in the 
shell, and thus must be ~ 10 km.s ^ Since the shell is a few par- 
sees wide, the shell matter needs to be renewed over a timescale 
of a few 10^ yrby on-going photo-evaporation. 

In the Eagle Nebula, the pressure from stellar winds is too 
low to account for the shell colors with collisional excitation (see 
section |5] for details). The mechanical power from the winds is 
also too small to contribute to the I R luminosity from the s hell. 
For a wind velocity of 2500 km.s ' jKudritzki & Pulsl2000l) . the 
mechanical energy injection is ~ 2500 Lq, a factor 20 smaller 
than the shell luminosity ~ 5 x 10^ Lq as estimated from the 
shell brightne ss Bjr and its an gular diameter (14'). Unlike what 
Everett & Churchwelll (1201 Oh advocated for N49, in Ml 6 the 
shell IR emission cannot be powered by the stellar winds, and 
does not represent a major cooling channel that impacts the dy- 
namical evolution of a wind-blown shell. 

The shell must originates from the only available source 
of dust, i.e., evaporating dense gas condensations within the 
ionization boundary of the Nebula. The difficulty in being 
certain that this is the right interpretation comes from inter- 
stellar dust (see sections|4]and|5](. Indeed, our dust modeUng 
in section |4] shows that the shell SED cannot be fit with the 
standard intersteUar dust size distribution. The fits shown in 
Figure [T4(d)| and [TSl iUustrate the uncertainty of the model- 
ing. It is beyond the scope of this paper to explore in a sys- 
tematic way the full range of possible solutions, but we are 
confident that any fit will involve shattering of dust grains to 
nanometric sizes. 

As a consequence of such an interpretation for the Eagle 
Nebula shell, we conclude that massiv e star forming r egions 
have a major impact on carbon dust. iGalliano et al.l (l2003l) 
reached a similar conclusion in their modeling of the infrared 
SED of the dwarf, star forming, galaxy NGC 1569. Observations 
of the ionized gas kinematics do provide evidence for supersonic 
velocitie s in the immediate environm ent of pillars in star forming 
regions (iWestmoquette et ani2009h . Hence, the grinding of the 
carbon dust could be the result of grain shattering in grain-grain 
collisions within shocks driven by the dynamical interaction be- 
tween the stellar winds and the shell. Theoretical modeling of 
the dust dynamics in s hocks suggest that t his is a plausible hy- 
pothesis (|Jonesll2004l) . lGuillet et alJ (l2009h have quantified dust 
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processing by the passage of J-shocks of a few 10 km.s They 
find that the mass fraction in the largest grains is reduced to the 
profit of the smallest, as a result of grain shattering and dust va- 
porization. 

6.2. A supernova remnant 

Alternatively, we keep the usual distribution of dust grain 
sizes, but look for another source of pressure: a supernova 
remnant. This is not unexpected for a 3-Myr old nebula with 
very massive stars (M* ~ 8OM0 Hillenbrand et al. 1993). If 
so, we would be witnessing a specific time in the evolution 
of the remnant where the plasma pressure and temperature 
would be such that the remnant cools through dust emission. 
This scenario relates directly to the fit of the shell SED quan- 
tified in section m 

The infrared dust emission from fast shocks driven by su- 
pernovae has been qua ntified in several theoretical papers (e.g. 
iDraind [T98TI: iDwek et al. 1996). Overall, dust is found to be 
a significant but not dominant coolant of shocked plasma due 
to dust destruction. For a dust to hydrogen mass ratio of 1% 
and a Solar metallicity, dust cooling is larger than atomic cool- 
ing for temperatures > 5 x 10^ K, but, for temperatures T 
larger than ~ 10^ K, the dust destruction tim escale by sput- 
tering is s maller than the gas cooling time (Smith et"al]|1996l: 
iGuillard et al. 2009). This framework has been used to interpret 
observations of young remnants starting from th e first infrare d 
detections of supernovae with the IRAS survey (lDweklll987h . 
We propose here a distinct idea, where the shell infrared emis- 
sion seen towards the Eagle Nebula would be related to the late 
evolution of a remnant. 

For the model shown in Fig. [18] 1/3 of the shell infrared 
emission is powered by grain collisions with electrons and con- 
tributes to the plasma cooling. The remaining 2/3 is provided 
by radiative heating of the dust. Assuming that the dust in- 
frared emission is the dominant gas cooling channel, the iso- 
baric cooling time of the infrared emitting plasma is tcooi - 
I X 2.3 X kTel(T X nip x xj) where F is the collisional heat- 
ing rate per unit dust mass, nip the proton mass and Xd the dust 
to hydrogen mass ratio. With the F value derived from the fit 
in Fig. [HI we find W = 1500 x (xd/O.Ol)"' yr. The dust- 
to-hydrogen mass ratio Xd is not constrained by the modeling. 
This factor may well be smaller than the reference value of 
1% due to dust destruction by sputtering. The SED fit also al- 
lows us to estimate the plasma column density and thereby the 
internal energy U of the infrared emitting plasma. The model 
gives Nh - 10'^ x (xrf/0.01)"' H.cm"^. From there we find 
U - 2 X 10"*** X (xd/O.Oiy^ erg. This value is a small fraction 
of the expansion energy associated with a typical supernova ex- 
plosion (~ 10^' erg). Within our remnant hypothesis, this large 
difference indicates that the cooling time is short and that only 
a small fraction of the shocked plasma is contributing to the in- 
frared emission. One possibility to account for this fact would be 
that we are observing the late evolution of the remnant when the 
low density hot plasma heated to high temperatures early in the 
expansion of the remnant is cooling through turbulent mixing 
with photo-ionized gas (Begelman & Fabian 1990). This plasma 
would have a long intrinsic cooling timescale, because its dust 
would have been destroyed early in the evolution of the remnant. 
For a pressure ofp/A: = 5xlO^ K.cm^^*, the cooling timescale 
through atomic processes of a dust-free plasma at a temperature 
of 10^ K is 2 X 10^ yr 

This interpretation will need to be tested against addi- 
tional observations. The absence of bright difi'use emission 



in the Chandra X-ray images dLinskv et alJl2007h can possi- 
bly be accounted for. For instance, the hot plasma may be 
too tenuous to be seen in emission, while the X-ray emis- 
sion from the turbulent mixing layers would be soft and thus 
heavily attenuated by foreground gas. We re-analyzed the 
Chandra ACIS-I observations of M16 (Linskyetal. 200^ 
to search for a faint background emission. After removal of 
point sources, we do find residual X-ray emission over the 
SW section of the mid-IR shell where the foreground extinc- 
tion is the lowest. The emission spectrum fit gives kT in the 
range 0.6 - 2 keV and a foreground column density within 
2.45.4 X 10^^ H.cm^. The absorption corrected X-ray bright- 
ness is 1.3 X lO"* erg.s'.cm^.sr'. If this emission arises from 
the mid-IR shell (i.e. from a sightline length ~ 10 pc), we de- 
rive a plasma pressure p/k ~ 10** K.cm^. This result does not 
allow us to conclude that the X-ray emission arises from a 
supernova remnant, but, if it does, the X-ray emission is con- 
sistent with the dust being coUisionally excited in a high pres- 
sure plasma. In this case, if the X-ray emission fills the mid- 
IR cavity, the shell X-ray luminosity would be ~ 10^^ erg.s'. 
This is on the low side for an SNR: for comparison, the W28 
SNR, which is interacting with a molecular cloud, has a to- 
tal X -ray luminosity ~ 6 x lO'^^erg s ' (Rho & Borkowskil 
I2OO2I) . However, our value of Lx for the putative M16 SNR is 
a lower limit, since it does not take into account the soft X- 
ray emission from cooler gas that is more heavily absorbed. 
Further X-ray observations are planned to clarify this point. 
MIR spectroscopic maps of M16 with Spitzer, covering a 
wide range of emission features and ionization energies, will 
provide an additional test to be investigated. 

7. Conclusions 

- We present new IR images of the Eagle Nebula from the 
MIPSGAL survey that reveal the well-known illuminated 
clouds of dust and gas. The MIPS 24 /im observations shows 
the same inner shell-like feature as mid-infrared observations 
from ISO or MSX. It is significantly brighter than the PDRs. 
Relative to these previous observations, the MIPSGAL sur- 
vey has the advantage to also probe the far infrared emission 
of the dust. The structure of the nebula as seen in the MIPS 
70 yum observations is close to that of the shorter wavelengths 
as seen in the GLIMPSE survey (from 3 to 8 yum): the cloud 
surface is significantly brighter than the inner shell. 

- Thanks to the MIPS 24 and MIPS 70 yum observations, we 
are able to give constraints on the temperature of the grains 
emitting in the FIR range and the required interstellar radia- 
tion field intensity to heat them up to these temperatures with 
our dust model. The dust temperature varies from ~35 K in 
the PDRs to ~70 K in the shell. The required intensity of the 
ISRF within the PDRs is about an order of magnitude lower 
than that provided by the star cluster NGC66 1 1 . The shell of 
hot dust, however, requires an ISRF intensity about a factor 
of a few higher than that provided by the cluster 

- Combining all the IR observations at our disposal into SEDs 
that sample the whole nebula with our dust model, we fit 
the observations to constrain both the radiation field inten- 
sity and the dust size distribution. In the PDRs, we confirm 
the required ISRF intensity is about a few tenth of that pro- 
vided by NGC661 1. The dust size distribution is dominated 
by BGs even though all the dust components are present with 
abundance a factor of a few, at most, different from those 
of the DHGL medium. In the shell, we also confirm the re- 
quired ISRF intensity is a factor of a few larger than that 
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of NGC661 1. The PAHs are absent and the VSGs are more 
abundant, up to a factor 10, than in the DHGL medium. 

- Extinction and the dispersion of the stars across the neb- 
ula can account for the lower ISRF intensity required for 
the PDRs. On the contrary, an additional source of heating 
is required for the shell. Neither the spatial variations of the 
ISRF intensity nor the Lyman alpha photons contribution can 
account for the discrepancy between required and provided 
UV heating of the dust. Exact positions of the stars reveal 
that Pilbratt's blob is only 0.25 pc from an 08.5V star and 
may thus be a bow shock. 

- We then invoke gas-grain collisions as an extra source of 
heating. Our modeling leads to a fit of the shell SED that 
requires a pressure of a few 10^ K.cm"^. Such a pressure is 
at least a factor of a few larger than that inferred either from 
optical observations at the end of the photo-evaporation flow 
arising from Pillar or from Pilbratt's blob bow shock nature. 

- We finally discuss two interpretations of the mid-IR shell 
in the general context of a massive star forming region. 
In a first scenario, we propose that the shell is wind 
blown by the stars. We find that the star cluster does not 
provide enough mechanical energy via stellar winds to 
power the shell emission. Therefore the shell is explained 
by a modified dust grain size distribution (large carbon 
grains shattered to nanometric sizes) with heating only 
due to UV emission. The implication is then that massive 
star forming regions like Ml 6 have a major impact on 
their dust size distribution: this can be checked on other 
similar regions. Alternatively, we propose a second sce- 
nario, in which the shell is heated by the hidden remnant 
of a supernova from a very massive progenitor, and for 
which the dust provides a fast cooUng. The impUcation 
is then that our observations occur during a short-Uved, 
late stage of evolution of the remnant: this can be checked 
with new X-ray observations. 

The Eagle Nebula IR emission morphology is similar to 
that of many other star forming regions observed within the 
GLIMPSE and MIPSGAL surveys ( Churchwell et al. 2006; 
ICarev et al.l2009h . For the first time, it is quantitatively discussed 
in terms of dust modeling. The work we present would need to 
be extended to other SFRs with IR morphology similar to that 
of Ml 6 to ascertain whether the interpretation would be chal- 
lenged by the same problem in accounting for the dust tem- 
perature. Moreover, future analysis of additional observations 
(mid-to-far IR spectral mapping from Spitzer/IRS and MIPS- 
SED, near-IR narrow band imaging from CFHT/WIRCam) of 
the Eagle Nebula will provide us with more constraints on the 
physical conditions and dust properties in M16's inner shell. 

This work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract 
with NASA. Support for this work was provided by NASA 
through an award issued by JPL/Caltech. 
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